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I . Introduction 


Recent studies of the radio emission spectrum of Venus 
[Steffes, 1986 and Steffes gt al. . 1990] have suggested that 
variations (both temporal and spatial) in the abundances of sulfur- 
bearing gases may be occurring below the Venus cloud layers. One 
technique for monitoring the abundances of sulfur-bearing gases 
(especially gaseous H 2 S0 4 ) is to measure the 13 cm wavelength 
opacity using Pioneer-Venus Orbiter Radio Occultation studies 
( PVORO ) [Steffes, 1985]. 

In order to characterize possible variations in the abundance 
and distribution of subcloud sulfuric acid vapor, we have processed 
13 cm radio occultation signals from 23 orbits that occurred in 
late 1986 and 1987 (Radio Occultation Season #10) and 7 orbits that 
occured in 1979 (Radio Occultation Season #1) . During Season #10, 
the a tmosphere was probed to altitudes as low as 39.4 km (above a 
mean radius , of 6051.92 km) at latitudes in the northern hemisphere 
from 11 to 89°N. During the orbits from Season #1, the “atmosphere 
was probed to altitudes as low as 36.7 km. The data were inverted 
to produce 13 cm absorptivity profiles. Finally, pressure and 
temperature profiles obtained with the Pioneer-Venus night probe 
and the northern probe [Seiff et al. , 1980] were used along with 
the absorptivity profiles to infer upper limits for vertical 
profiles of the abundance of gaseous H 2 S0 4 . In addition to 
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inverting the data, error bars have been placed on the absorptivity 
profiles and H 2 S0 4 abundance profiles using the standard propagation 
of errors [Brandt, 1963]. These error bars were developed by 
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considering the effects of statistical errors only. The profiles 
show a distinct pattern with regard to latitude which is consistent 
with latitude variations observed in data obtained during 
occultation seasons #1 and #2 (late 1978 and early 1979) [Cimino, 
1982]. However, when compared with the earlier data, the recent 
occultation studies suggest that the amount of sulfuric acid vapor 
occurring at and below the main cloud layer may have decreased 
during the time interval between early 1979 and late 1986. 

II. Processing of Radio Occultation Data 

Radio occultation studies using the Pioneer-Venus Orbiter 
allow us greater flexibility in studying the vertical structure of 
the neutral atmosphere of Venus than either in-situ probes or radio 
astronomical observations. The atmosphere can be probed at a wide 
range of latitudes due to the geometry of the spacecraft's orbit 
and its changing orientation with respect to Earth throughout an 
occultation season. In addition, the longevity of the orbiter has 
made it possible to conduct long term studies of the microwave 
properties of the atmosphere of Venus with greater vertical 
resolution than radio astronomical observations. 

Complete descriptions of the theoretical bases for this 
technique are given in Fjeldbo and Eshleman (1965) , Fjeldbo al.. 
(1971) and Eshleman (1973). For the Pioneer-Venus application, 
consider the geometry illustrated by Figure 1. In this experiment, 
the spacecraft modulation is switched off and a continuous wave 
(CW) signal is transmitted during the occultation event. For a 
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successful experiment, the radiated carrier signal must have a 
stable amplitude and frequency. To ensure that the frequency is 
stable, a reference signal is uplinked to the spacecraft from the 
Earth station. The spacecraft locks its oscillator to this signal 
and downlinks a noninteger harmonic of the received frequency. The 
power output of the spacecraft is stabilized internally. The 
occultation experiment begins when the path of the signal 
transmitted by the spacecraft tangentially grazes the edge of the 
atmosphere of Venus. As the spacecraft's trajectory carries it 
toward the limb of the planet, the path of the signal slices deeper 
into the atmosphere, traversing a longer pathlength through the 
limb. Refraction in the atmosphere causes the signal to bend 
around the limb of the planet. This increases the Doppler shift 
since a larger component of the velocity of the spacecraft is now 
parallel to the ray path back to Earth. It is important to note 
that this also affects the uplinked reference signal from Earth. 
Hence the Doppler shift of the signal received at Earth includes 
the shift incurred by the uplink to the spacecraft as well as by 
the downlink. 

The power of the received signal is attenuated by two 
atmospheric effects and by mispointing of the spacecraft antenna. 
The curvature of the atmosphere and the changing refractive index 
profile cause spreading of the transmitted beam, an effect known as 
refractive defocusing. Secondly, absorption and scattering in the 
atmosphere extract energy from the signal, causing an additional 
drop in received power. In addition, as the ray is bent by 
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refraction in the atmosphere, it is also bent away from the main 
axis of the spacecraft antenna, causing a drop in power due to the 
smaller off-axis gain of the antenna. As the experiment 
progresses , the Doppler shift increases and the amplitude drops 
until the spacecraft can no longer maintain lock with the uplinked 
frequency reference. 

For Figure 1, we assume that for each instant in time, the 
illustrated distances R e , D and y (and hence, R es , R s and the angle 
y are known from interpolations of spacecraft tracking data, from 
Earth and planetary ephemerides and from the knowledge of the 
tracking station position relative to the center of the Earth. Thus 
the other quantities follow as: 

a = Dsin (6 - P) + ycos(6-P) 

6 = (6-P) + sin' 1 (a/R e ) 

( 1 ) 

R x = (Rl - a 2 ) 1/2 
r 2 = (j?J - a 2 ) 1/2 


[Eshleman et al. , 1980] . 

The impact parameter a and the bending angle 6 can be used to 
calculate the vertical profile of refractivity if we assume a 
perfectly stratified atmosphere with no horizontal variations other 
than planetary curvature in the region probed by the beam [Fjeldbo 
et al . , 1971] : 


In (n(r) ) 


1_ f m 6 (a) da 

X Ja(r) (j,) 


( 2 ) 
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The lowest depth, r, probed for each ray is related to the ray 
impact parameter by a = rvr. Thus, the natural log of the index of 
refraction n is the inverse Abel transform of the integral of the 
bending angle <5. Once n(r) is determined it is possible to obtain 
pressure and temperature profiles, assuming hydrostatic equilibrium 
and that the major constituents are well known and well mixed [Lipa 
and Tyler, 1979]. 

In order to invert the amplitude data to obtain absorptivity 
profiles it is necessary to first remove the effects of refractive 
defocusing and mispointing of the spacecraft antenna. The 
refractive defocusing can be computed using geometric optics and is 
then subtracted from the signal attenuation profile. Let I be the 
received signal intensity normalized by the free space intensity at 
the spacecraft receiver input. Then the reduction in intensity due 
to refractive defocusing, I' 1 , is given by 


y x / cos (6-P) dy v 
R es sin P R es d P 


(3) 


where the first bracketed term represents the focusing effect of 
the curved limb normal to the plane of the figure and the second 
term is the refractive defocusing in the plane of the figure 
[Eshleman £t &1. , 1980]. For an Earth-Venus distance much larger 
than the spacecraft-Venus distance, equation (3) reduces to 


I - 1 = 


a <i 


( 4 ) 


5 


[Eshleman £t al . , 1980] . 

To understand the importance of the antenna pointing issue, it 
must be noted that, as currently operated, the direction of the 
beam of the High Gain Antenna (HGA) aboard the Pioneer-Venus 
Orbiter is only adjusted once per orbit so that it is pointed 
directly to Earth at the periapsis of each orbit. Since the ray 
path back to Earth can be bent by as much as 3.5° from the true 
direction to Earth, significant signal reductions occur. These 
reductions are not due to atmospheric effects, but are due to the 
gain dependence of the orbiter antenna on incidence angle. During 
the Pioneer-Venus primary mission, slew commands were sent to the 
orbiter so that the antenna beam would continually track the 
virtual image of the Earth. Thus, little if any correction was 
needed for power variation due to S-band antenna mispointing. 
However, since such tracking was not used during Season #10, it is 
necessary to trace the ray path for each data point and determine 
the gain of the orbiter antenna at the ray's angle of departure. 

After correcting the amplitude data for the effects of 
refractive defocusing and for antenna mispointing, we are left with 
the excess attenuation profile r . This residual attenuation is the 
result of integrating the absorption coefficient a over the 
downlink ray path for each sample point. (This assumes that 
scattering is negligible compared to absorption at the transmitted 
frequency.) This can be written as an integral over radius using 
Bouger's rule [Lipa and Tyler, 1979]: 
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a (r) dr 


( 5 ) 


T U 0 > 




a(-r 0 ) r 

r*ii(r) ' 
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where r is in dB and a is expressed in dB/km* This expression can 
be rewritten in the form of an Abel transform by making the 
approximation da « n-dr: 


t u 0 > 


O f « (r(a) ) 
Jr 0 n(r{a ) ) 


a da 

y/a 2 - a (r 0 ) 2 


( 6 ) 


The form of the Inverse Abel Transform that has been applied to the 
data is the following: 


£(£o) . dj 2 f~ t (a) a da 

2 a (r 0 ) dal J aUo) ^ a 2_ a ( Xo )2 


( 7 ) 


Thus, the absorption coefficient can be recovered as a function of 
height by applying an inverse Abel transform to the excess 
attenuation r and multiplying the result by the refractivity 
profile. 

Although this is a fairly complete description of the 
mathematical procedure for reducing radio occultation data, a 
complete analysis requires characterization of experimental 
uncertainties. Since the derived quantities a and n cannot be 
observed directly, neither can the uncertainties associated with 
each of them. The next section describes our method for placing 
error bars on the absorptivity profiles. 
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III. Error Analysis 


Statistical errors are inherent in any set of experimental 
observations. The significance of small scale structure in 
experimental data is dependent on the magnitude of these 
statistical errors. Thus, knowledge of the variance of the noise 
fluctuations riding on the true data is necessary to evaluate the 
validity of theoretical predictions or to compare independent sets 
of observations. In the case of indirect observations, where the 
desired vector quantity, y, is obtained by transforming the 
observed vector of quantities, x, it is not possible to measure the 
noise power directly. If y is related to x by a linear transform 
T , then the covariance matrix of y is given by 

yx 

C - T C T T (8) 

c y 1 yx 1 yx 


where 

C x = < (x-x) -(x-£) T > (9) 

[Brandt, 1963]. Here, x is the vector of true values, and £ is the 
measured, or average value and <•> is the expectation operator. If 
the errors are small and the transform T^ is nonlinear, then C y can 
be estimated by expanding T yx in a Taylor series 

y = t 0 + T x (x-£) + T 2 (x-£) 2 + ... (1°) 


where 
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d Yi | 

d Yi | 

d Yi 

1 

dx x ]x =* 

dx l*=* - 

dx n 

lx=je 

dy 2 1 

d y 2 1 

d y 2 

1 

dx x lx ‘* 

dx 2 ' x ‘* - 

dx n 


d Vm I 

d y* | ... 


1 

dx x lxm * 

dx 2 1 x “ * 

dx n 

lx-je 


Now the estimate for the C y is given by 

c y = r 1 -c x T 1 r 


(ii) 


( 12 ) 


This is the result of neglecting all terms higher than first order 
and is known as the standard propagation of errors [Brandt, 1963]. 

To see how this method may be applied to a multistep procedure 
such as the reduction of radio occultation data, consider the 
vector quantity z = T 2 y. Now, C z is given by 

c z = t 2 c y t 2 


= t 2 {t x c x t*) t 2 t 

= (T 2 T X ) C X (T 2 T X ) T 


( 13 ) 


= r,c v 7Y 


where T 3 = T 2 T 1 . This forms the basis for placing error bars on the 
absorptivity profiles discussed in Section IV. Note that this 
procedure can be carried out with fewer numerical computations if 
the sequence of transformation matrices are multiplied in the 
appropriate order before the covariance matrix for the desired 
quantity is determined. For example, if k transformations are used 
which can be represented as n x n matrices, then forming the 
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covariance matrix at each step (C i = ’*VC j1 , T i T ) requires 2n 2 
multiplies for a total of 2 Jen 2 operations. By contrast, forming 
the product T f = T k -T k1 * . . . • Ty T, requires kn 2 multiplies and forming 
the covariance matrix C = TycyT, 1 requires an additional 2n 2 
multiplies, resulting in only (k+2)n 2 total operations. 

Figure 2 illustrates the propagation of errors through the 
various steps of inverting the frequency and amplitude data to 
obtain the absorptivity profile. In the figure, each 

transformation of variables is represented by a labeled box. The 
random variables which flow from each box are related to the random 
quantities that flow into that box by that transform. For example, 
the ray path parameters are determined from the frequency data and 
trajectory information. Thus, A a and A S are dependent on the 
uncertainties in the frequency data, Af. In the first step, errors 
are estimated for power and frequency. Lipa and Tyler (1979) 
developed a model for estimating these errors based on their method 
for estimating the amplitude and frequency of the sampled signal 
received from the spacecraft. In the case of this study, the 
available data consisted only of the excess Doppler shift (Hz) and 
relative power (dB) as time series. Because the original sampled 
signal is needed to estimate the errors in frequency and power by 
Lipa and Tyler's method, an alternative method has been developed. 

The approach used is to fit a line to the set of points 
{ f (t n . 4 ) f f (t n . 3 ) , . . . , f (t n _.,) , f (t n ) , f (t m1 ) , • • • / i (t^j) , f (t^) } for each 
point t n in time, and use the mean square error as the estimate for 
Af n , the variance for the frequency at point n. The same procedure 


10 


is used to estimate the variance at each point of the power 
profile. Since the noise process is very nearly white Gaussian 
noise and is stationary over short intervals [Lipa and Tyler, 
1979], this provides a simple method for estimating the errors Af 
and Ap. We are assuming (as did Lipa and Tyler) that the random 
processes Af and Ap are independent and uncorrelated so that 
<p M - f n > = 0 for all m and n, and <p m ‘ p n > = <f m ' f n > = 0 for m * n. 
Thus the covariance matrices C f and C p are diagonal and the 
covariance matrix C fp is identically zero. 

The next step involves determining the uncertainties in the 
ray path parameters, which are derived from the frequency data and 
the spacecraft trajectory. The equations used to determine these 
parameters are linearized and the standard propagation of errors 
used so that 

= TifCfTif 

C a = T af C f Tj f (14) 
C 8a = T hf C f T? f = C a l 

Since C f is diagonal, C 5 , C 5a and C a are also diagonal at this point 
of the processing. This assumes that the uncertainties in the 
trajectory data are negligible compared to the uncertainties 
arising from the noise on the received signal. 

Before the index of refraction is calculated, all profiles, 
such as a, 6, and p are "averaged and decimated" by partitioning 
the points of each profile into bins based on ray asymptote height. 
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The size of the bins is specified by the user and, in roost cases in 
this study, is 0.5 km. The points in each bin are replaced by a 
single point which is the average value of all the points in that 
bin. This is a useful technique since the derivative d<5/da 
increases as the occultation progresses, causing the points to 
become denser in altitude with respect to time. The result is a 
set of vectors which have been decreased from a typical length of 
over 400 to less than 100 points. An additional benefit is that 
the variance of each new point is roughly a factor of 1/n times the 
variances of the n points in that bin. Also, the points at the 
bottom of each profile, which are more noisy than upper points, 
have experienced more averaging, and hence, greater reduction in 
noise, due to the fact that the lower bins tend to have more points 
in them. 

Because of this significant reduction in the lengths of all 
the vectors used in the remaining data reduction, it is convenient 
to trace the errors as they propagate through the data processing 
by considering the transformations of the vectors <5, a and p. 
Thus, all references to the covariance matrices associated with 
these vectors, C s , C a , or C p , or the uncertainties A<5, Aa or Ap are 
with regard to this stage of the processing. 

The uncertainties associated with the index of refraction is 
determined as per Lipa and Tyler (1979). The exponential function 
imbedded in equation (2) is approximated by linearizing about zero, 
the integral is approximated by a summation using the trapezoidal 
rule and the result is linearized with respect to the variables <S 
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and a so that 


An = r n4 A6 + T na Aa 


(15) 


where is the linear transform relating 6 to n (to the first 
order) , and T M is the transform relating a to n. 

From the fact a = n-r, it follows that 


Ar = T za Aa + T rn An 

= (Tra + Trn T n a )ba+T in T n6 Ab 


(16) 


where T ra is a diagonal matrix with diagonal entries (T ra ) i f = n/ 1 and 
T m is a diagonal matrix with diagonal entries (T rn ) H = -a./n 2 . 

After the calculation of the index of refraction, all vectors 
are averaged and decimated again, only this time with respect to 
the radius profile instead of the ray asymptote. Then, if desired, 
the data is smoothed by a 5 point running average. Let T ave denote 
the matrix representing the process of averaging and decimating the 
vectors, and T s denote the matrix representing the smoothing 
process . 

Now the excess attenuation is determined by performing the 
antenna correction and removing the refractive defocusing. We have 
determined that the uncertainties in performing the antenna 
correction are small compared to the other sources of error (by 3 
orders of magnitude) and have chosen to neglect them. We have 

Ax = r Tp Ap + T xa T s T av0 Aa + T xh T s T ave A 5 (17) 
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where T and T . are due to refractive defocusing. 

T0 TO 

The inverse Abel transform of the excess attenuation is 
computed in two steps. Examining equation (7) reveals that the 
integral corresponds to the forward Abel transform of r . The 
result of this transform is arbitrarily named F and has the 
associated uncertainty 

A F = TVjAt + T Fa Aa 

« T ft T tp Ap + T^T xh T a T av ^b 

( 18 ) 

+ ( Tp, r Ta T s T ave + T Fa ) A a 
= T'fp A p + T'pi Ad + l> a Aa 


The derivative of F is taken with respect to a and scaled by the 
factor -n/(27ra) to yield the absorptivity profile a. The 
uncertainty in a is given by 

Aa = T aF LF + T aa ^a + r #n An 

= + 7> 6 A6 + T' Fa Aa) 

+ T aa A a + T an {T nb A6 + r na Aa) 

(19) 

= TaF T ?p A P + ( T *F T Fl + Tan T n 6 > A6 
+ (T uF T Fa + T aa + T an r M ) Aa 
= r op Ap + r ## A8 + r #a Aa 

Finally, it is possible to compute the covariance matrix of a: 
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<Aa ’Aa T > 


C a = 

- T. p C p T.Z * T. s C,tJ * T..C.T.J ( 20 ) 

+ ^6 a ^aa + ^«a 

Since C tf8 = C^ T , we need only compute one of the final two terms in 
the s um above. The error bars displayed on the figures discussed 
in the next section are the square roots of the diagonal elements 
of the matrix C a . 


IV. Discussion of Results 

Table I lists the 23 data sets from Season #10 (July 23, 1986 
through January 7,1987) and the 7 data sets from Season #1 
(December 28, 1978 through January 26, 1979) which have been 
reduced by the method described in section II. The latitude and 
solar zenith angles for the deepest point probed in each experiment 
are given, as well as the maximum absorptivity measured before loss 
of lock (in the case of the entry occultations) . Most of the 
occultations were entries and are denoted by the suffix "N". Orbit 
2787X is the only data set analyzed which resulted from an exit 
occultation. We found that the data sets from most exit 

occultation experiments had to be discarded due to the instability 
of the spacecraft's onboard oscillator. Although the frequency 
data from 2787X appears to be relatively stable (in that the data 
could be processed) , the results for this experiment are not deemed 
to be as reliable as the data from the entry occultations. 

Figures 3 through 5 show the absorptivity measured for entry 
occultation experiments conducted for orbits 2787, 2801 and 2928, 
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respectively. Figures 6 through 10 show averages of the 13 cm 
absorptivity profiles we have currently obtained for 5 different 
latitudinal areas, based on the 23 profiles currently available. 
Figure 6 is the result of averaging profiles 2787N and 2792N. 
Figure 7 is the result of combining 2798N and 2801N. The profile 
in Figure 8 is the result of averaging 2814N, 2815N, 2819N, 2859N, 
2860N, 2862N and 2952N. Profiles 2844N, 2845N, 2850N, 2851N and 
2853N were combined to produce Figure 9. Figure 10 is the result 
of averaging 2921N, 2923N, 2928N and 2930N. Profiles 2939N and 
2955N were not included in the averaged profiles of Figures 9 and 
8, respectively, due to the relatively high altitude at which loss 
of signal occurred in these experiments. Profile 2787N was not 
included due to the unknown systematic errors resulting from the 
unstable spacecraft frequency for exit occultations. 

When the occultation profiles for the different areas are 
compared, it appears that the level of significant absorptivity for 
the various latitudinal regions decreases with increasing latitude. 
A similar effect had been observed with the 1979 data, although the 
difference between the levels for the polar occultations and the 
collar occultations (latitudes in the 70° s) was severe. Since the 
13 cm opacity is principally due to gaseous H 2 S0 4 [Steffes, 1985], 
this suggests a variation in the altitude of the cloud base with 
latitude (i.e. the altitude below which the absorption occurs). 
Although, the error bars are quite large (there is no statistically 
significant absorptivity in many of the profiles) , we view the 
trends in the latitudinal variations as being valid. This is 
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especially noticeable in Figure 11, where the average absorptivity 
profile for the equatorial zone is compared with that for polar 
latitudes . 

Moreover, when these absorption profiles (and the values for 
peak absorptivities below the clouds, as shown in Table I) are 
compared with those obtained during the initial Pioneer-Venus radio 
occultations (1978-79) [Cimino, 1982], it appears that the average 
13 cm opacity has dropped, implying a reduced abundance of sulfuric 
acid vapor. While this is suggestive of long-term temporal 
variation in the subcloud abundance of gaseous H 2 S0 4 , it cannot be 
considered proof until the uncertainties in the earlier 
absorptivity profiles are properly characterized. We have newly 
processed 7 orbits from Radio Occultation Season #1 (including 
error bars) . An example of these profiles is shown in Figure 12 
(Orbit 40N) . The larger significant absorptivity is noteworthy, 
especially when compared with the average for the 70° latitude 
range obtained during Season #10 (Figure 9) . An average of orbits 
from Season #1 at latitudes above 70°N (24N, 40N and 44N) is given 
in Figure 13, along with the averaged profile for 70° obtained for 
Season #10. This figure suggests graphically that the abundance of 
H 2 S0 4 decreased from 1979 to 1986, although the error bars for both 
profiles do overlap somewhat. 

Figure 14 through 36 display the abundance of sulfuric acid 
vapor as a function of radius inferred from the derived 
absorptivity profiles from all 23 orbits analyzed from Season #10 
and the saturation vapor pressure using results from a recent 
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laboratory measurement by Fahd and Steffes (1990) . For Figures 14 
through 18 , the pressure and temperature measurements from the 
Pioneer-Venus night probe were used to derive the abundance 
profiles. The pressure and temperature measurements from the 
Pioneer-Venus north probe were used for Figures 19 through 36 
[Seiff, Si . , 1980]. These figures show the extent to which 
interpretation of the amplitude data from occultation experiments 
can be carried, that is, to determine upper limits for the 
concentration of opaque constituents in planetary atmospheres. 

V. Summary of Activities 

In the 27 months of Grant NAG 2-515 (July 1, 1988 through 
September 30, 1990) we have worked on reducing 13 cm radio 
occultation data to obtain 13 cm absorptivity profiles and the 
resulting abundance profiles for gaseous H 2 S0 4 . A summary of our 
activities is given in Table II. During this period of time all 
previously existing software was corrected and modified to apply 
the standard propagation of errors to each step in the processing 
of the data to derive error bars for the absorptivity profiles and 
H 2 S0 4 abundance profiles. This activity has been conducted using 
the JPL-Radio Occultation Data Analysis Network (JPL-RODAN) , which 
was accessed remotely from Georgia Tech both by telephone modem and 
by INTERNET/ ARPANET. We worked with JPL to improve the function of 
its INTERNET interface. This allowed us to off-load some data 
processing and data analysis to Georgia Tech machines. Graduate 
student Jon M. Jenkins and Professor Steffes visited JPL 5 times 
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during the course of the project to load the data sets from 
magnetic tape and process them with the new error estimating 
software. Also, we attended 5 Pioneer-Venus Science Steering Group 
(PV-SSG) meetings and have worked effectively with the Radio 
Science Team Leader (Dr. A. J. Kliore) to evaluate the data. 

In May 1990, we submitted a paper to the journal I<?apu§ 
entitled, "Results for 13 cm Absorptivity and H 2 S0 4 Abundance 
Profiles from the Season 10 (1986) Pioneer~Venus Radio Occultation 
Experiment" [Jenkins and Steffes, 1990]. In this paper, we 
described our work in developing 13 cm absorptivity profiles and 
their accompanying error bars for the orbits listed in Table I. 
(Note that the 13 cm abundance profiles are directly related to the 
abundances of gaseous H 2 S0 4 , as discussed previously.) 


V. Conclusions and Suggestions for Future Work 

As noted, the error bars placed on the absorptivity profiles 
derived using the current method are quite large. This may be due 
in part to overestimating the uncertainties in the frequency and 
power profiles at the beginning of the process. However, we feel 
that it is better to overestimate the errors rather than to 
underestimate them. It is clear that the final step of 
transforming the excess attenuation into absorptivity via the 
inverse Abel transform introduces the greatest statistical errors 
into the measurements. This is due to the fact that the inverse 
Abel transform corresponds to half order differentiation 
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[Bracewell , 1977]. Thus, it may be possible to reduce the 
magnitude of the estimated errors and retain small-scale structure 
by employing an alternative method for evaluating this recalcitrant 
transform. 

In order to accomplish this task, software is being developed 
to simulate the occultation experiment. Given a model atmosphere 
and spacecraft trajectory, the software will generate time series 
for the received power and frequency. Noise will be added to these 
quantities and the results fed to the data inversion software. By 
inserting various methods for computing the inverse Abel transform 
into the inversion software, their performances can be objectively 
compared to select the method which is most stable and accurate 
under low signal-to-noise conditions. This simulation software 
will not only allow us to determine the method which best handles 
noise, but will also allow us to locate and evaluate some sources 
of systematic errors within the data inversion software. Graduate 
Assistant Jon M. Jenkins will continue to pursue this in his 
dissertation research. His Ph.D. Thesis (expected 6/91) will 
include both the Venus atmospheric results described here and the 
new techniques for processing radio occultation data. 

The Pioneer-Venus Orbiter Radio Occultation experiments 
represent a potentially powerful tool for examining the middle 
atmosphere of Venus. Our results suggest the presence of 
latitudinally-dependent variations in the abundance of sulfuric 
acid vapor and a global decrease in the overall abundance of H 2 S0 4 
below the main cloud layer. In the future, we hope to be able to 
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verify the statistical validity of these possibilities by modifying 
the techniques applied in the data inversion software and 
reprocessing additional data from earlier occultation seasons. 
Likewise, we hope that radio occultation opportunities using the 
more powerful Magellan spacecraft can be used to probe to deeper 
altitudes with higher accuracies than have been previously 
achieved. 
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TABLE I. DATA SUMMARY 


Orbit* 

Latitude 

(o) 

S Z A 
(°) 

Depth** 

(km) 

Maximum Absorptivity 
Measured Before Loss 
of Signal (dB/km) 

Season 10: 






2787N 

11.1 

160.6 

47.0 

0.0095 

± 

0.0035 

2787X 

-37.7 

139.9 

44.4 

0.010 

± 

0.001 

2792N 

24.8 

151.7 

45.9 

0.010 

± 

0.0045 

2798N 

36.1 

142.1 

43.7 

0.0085 

+ 

0.009 

2801N 

40.8 

137.7 

42.0 

0.0031 

± 

0.0028 

2814N 

59.3 

118.6 

41.3 

0.0036 

+ 

0.0039 

2815N 

60.7 

117.2 

40.6 

0.0044 

± 

0.0045 

2819N 

65.6 

111.6 

40.2 

0.0030 

± 

0.0049 

2844N 

79.6 

79.1 

40.9 

0.0054 

± 

0.0051 

2845N 

79.3 

78.4 

41.7 

0.0027 

± 

0.0055 

2850N 

75.6 

73.5 

41.9 

0.004 

± 

0.005 

2851N 

74.4 

72.2 

41.5 

0.0034 

± 

0.0050 

2853N 

72.6 

70.3 

41.2 

0.005 

± 

0.005 

2859N 

67.3 

65.3 

40.7 

0.0034 

± 

0.0051 

2860N 

66.7 

64.8 

41.0 

0.0066 

± 

0.0027 

2862N 

65.9 

64.8 

40.7 

0.0037 

± 

0.0057 

2921N 

87.2 

89.8 

40.5 

0.0064 

± 

0.0050 

2923N 

88.0 

90.9 

42.3 

0.001 

± 

0.007 

2928N 

88.6 

93.8 

39.4 

0.024 

± 

0.010 

293 ON 

87.5 

94.6 

40.6 

0.000 

± 

0.002 

2939N 

77.8 

104.0 

43.5 

0.0082 

± 

0.0019 

2952N 

60.4 

121.9 

42.0 

0.0042 

± 

0.0041 

2955N 

56.4 

126.0 

43.8 

0.0031 

± 

0.0029 

Season 1 

: 






24N 

87.0 

95.0 

39.1 

0.010 

± 

0.001 

4 ON 

78.0 

104 . 6 

40.9 

0.014 

± 

0.0014 

44N 

74.7 

107.8 

39.9 

0.013 

± 

0.007 

45N 

74.1 

108.4 

43.0 

0.0092 

± 

0.0019 

5 IN 

65.7 

116.6 

36.7 

0.011 

± 

0.007 

52N 

66.0 

116.2 

39.3 

0.0083 

± 

0.0017 

53N 

63.5 

118.7 

38.0 

0.0079 

± 

0.0007 


* Note that an N suffix indicates an entry occultation, and an X 
suffix indicates an exit occultation. 


** Deepest altitude probed (relative to a radius of 6051.92 km) 
before loss of signal. 
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TABLE II 


PIONEER— VENU S RADIO OCCULTATION (ORO) DATA REDUCTION: 
PROFILES OF 13 -CM ABSORPTIVITY 

FINAL REPORT 

1. TO DATE, 23 ORBITS HAVE BEEN PROCESSED FROM OCCULTATION SEASON 
10 (LATE 1986 AND EARLY 1987), INCLUDING ERROR BARS AND 7 
ORBITS HAVE BEEN PROCESSED FROM OCCULTATION SEASON 1 (EARLY 
1979) . 

2. PROBLEMS IN PREVIOUSLY EXISTING SOFTWARE HAVE BEEN DETECTED IN 
AREAS OF ANTENNA POINTING CORRECTION, SIGNAL GAIN 
NORMALIZATION (ALSO KNOWN AS AGC ACTION) , AND SIGNAL PATH 
DETERMINATION. CORRECTIONS HAVE BEEN DEVELOPED. 

3. PROFI LE S OF ABSORPTIVITY AND GASEOUS H 2 S0 4 ABUNDANCE AND THE 
ACCOMPANYING ERROR BARS HAVE BEEN DEVELOPED (2.29 GHz OR 13 
cm) . 

4. PROFI LE S GENERALLY SHOW A LOWER 13 cm ABSORPTIVITY THAN 
TYPICALLY MEASURED IN THE FIRST RADIO OCCULTATION SEASON 
(1979-1980) , HOWEVER, PROPER CHARACTERIZATION OF UNCERTAINTIES 
(I.E. ERROR BARS) OF EARLIER PROFILES IS NECESSARY TO CONFIRM 
ANY LONG TERM TEMPORAL VARIATION. 

5. DATA DOES NOT SUGGEST ANY DIURNAL VARIATIONS TO MATCH THOSE 
OBSERVED BY EMISION MEASUREMENTS AT 2 . 6 mm. 

6. LATITUDINAL VARIATIONS ARE APPARENT, I.E. ABSORPTIVITY BEGINS 
AT LOWER ALTITUDES IN POLAR REGIONS. 

7. ABUNDANCES OF GASEOUS H 2 S0 4 IN THE 5-20 PPM RANGE ARE COMMON 
BELOW THE CLOUD LAYERS. 

8. FUTURE WORK: 

INVESTIGATE ALTERNATIVE METHODS FOR COMPUTING THE INVERSE 
ABEL TRANSFORM SO AS TO REDUCE UNCERTAINTY IN DERIVED 
ABSORPTIVITY PROFILES. 

RADIO OCCULTATION STUDIES WITH MORE POWERFUL MAGELLAN 
SPACECRAFT. 
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Figure 1: Occultation Geometry 
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Figure 2: Flow chart for propagation of errors. 
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Figure 4: Atmospheric absorptivity profile measured for Orbit 

2801N which occurred on August 6, 1986 at 40.8°N. 
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Figure 6: Average 13 cm opacity profile obtained for the 

equatorial latitudes (11 to 25®N) . 
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Figure 10: Average 13 cm opacity profile obtained for the 

polar region (latitudes above 85°N) . 
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Figure 13: Average 13 cm opacity profiles for latitudes in the 

70° s from Season 10 (solid line) and latitudes above 
74°N from Season 1. 


38 



RROIUS . KM 







RADIUS . KM 


61151 


ORBIT 2787N 


6110 


6105 


610 0 



6095 | 

-20 


0 

H2S04 


_j i 

2 0 4 0 


RBUNDANCE , PPM 


6 0 


Figure 15: Sulfuric acid vapor abundance profile measured for 

orbit 2787N, which occurred on July 23, 1986 at 11.1°N 
(solid line) and the abundance expected from saturation 
vapor pressure (dotted line). 
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Figure 17: Sulfuric acid vapor abundance profile measured for 

orbit 2798N, which occurred on August 3, 1986 at 

36.1°N (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 
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Figure 18: Gaseous H 2 So 4 abundance as a function of radius 

measured for orbit 2801N, which occurred on August 6, 
1986 at 40. 8° N (solid line) and the abundance expected 
from saturation vapor pressure (dotted line) . 


43 



RADIUS » KM 


ORBIT 281 4N 


6115 


6110 


6105 


6100 


6095 


6090 



-40 


H2S04 ABUNDANCE 9 PPM 


Figure 19: Sulfuric acid vapor abundance profile measured for 

orbit 2814N, which occurred on August 19, 1986 at 

59.3°N (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 
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Gaseous H 2 So 4 abundance as a function of radius 
measured for orbit 2815N, which occurred on August 20, 
1986 at 60.7°N (solid line) and the abundance expected 
from saturation vapor pressure (dotted line) . 
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Figure 21: Sulfuric acid vapor abundance profile measured for 

orbit 2819N , which occurred on August 24, 1986 at 

65.6°n (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 
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Figure 22: 


Gaseous H 2 So 4 


abundance as 


function of radius 


measured for orbit 2844N, which occurred on September 
18, 1986 at 79 . 6°N (solid line) and the abundance 
expected from saturation vapor pressure (dotted line) . 
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Figure 23: Sulfuric acid vapor abundance profile measured for 

orbit 2845N, which occurred on September 19, 1986 at 
79. 3° N (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 
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Gaseous H 2 So 4 abundance as a function of radius 
measured for orbit 2850N, which occurred on September 
24, 1986 at 75.6°N (solid line) and the abundance 

expected from saturation vapor pressure (dotted line) . 
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Figure 30: Gaseous H 2 So 4 abundance as a function of radius 

measured for orbit 2921N, which occurred on December 
4, 1986 at 87.2°N (solid line) and the abundance 

expected from saturation vapor pressure (dotted line) . 
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Figure 31: 


Sulfuric acid vapor abundance profile measured for 
orbit 2923N, which occurred on December 6, 1986 at 

88.0°N (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 
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Figure 33: Sulfuric acid vapor abundance profile measured for 

orbit 2930N, which occurred on December 13, 1986 at 

87.5°N (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 
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Sulfuric acid vapor abundance profile measured for 
orbit 2939N , which occurred on December 22, 1986 at 
77.8°N (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 
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Gaseous H 2 So 4 abundance as a function of radius 
measured for orbit 2952N, which occurred on January 
4, 1986 at 60.4°N (solid line) and the abundance 

expected from saturation vapor pressure (dotted line) . 
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Figure 36: Sulfuric acid vapor abundance profile measured for 

orbit 2955N, which occurred on January 7, 1986 at 

56.4°N (solid line) and the abundance expected from 
saturation vapor pressure (dotted line) . 


61 



